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We investigate the gravitational implosion of magnetized matter by studying the inhomogeneous 
collapse of a weakly magnetized Tolman-Bondi spacetime. The role of the field is analyzed by looking 
at the convergence of neighboring particle worldlines. In particular, we identify the magnetically 
related stresses in the Raychaudhuri equation and use the Tolman-Bondi metric to evaluate their 
q , impact on the collapsing dust. We find that, despite the low energy level of the field, the Lorentz 

Pj^ ■ force dominates the advanced stages of the collapse, leading to a strongly anisotropic contraction. 

In addition, of all the magnetic stresses, those that resist the collapse are found to grow faster. 
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I. INTRODUCTION 

Magnetic fields are common features of almost all astrophysical environments and stellar magnetism is a long 
established and very active branch of astrophysics p}. Nevertheless, analytical studies of magnetic fields in strong 
gravity environments are less developed. Most of the available work addresses the possible gravitational effects on 
(j \ the Maxwell field (e.g. see and references therein) and relatively few look into the implications of magnetic fields 
O" 1 for gravitational collapse itself. Perhaps the most intriguing result so far has been obtained by Thorne in his analysis 
| | of Melvin's cylindrical magnetic universe There, by developing the concept of 'cylindrical energy', the author 
bi). reached the conclusion that "a strong magnetic field along the axis of symmetry may halt the collapse of a finite 
cylinder before the singularity is reached" [5[ . The possible support of the field against the gravitational collapse of 
massive bounded systems was also studied in [j| . That analysis led to solutions of Einstein-Maxwell equations with no 
singularities or event horizons, where the gravitational attraction is balanced solely by magnetic stresses. Other work, 
however, indicated that locally naked singularities may develop in some cases [6|. Studies of contracting charge dust 
has suggested that the fluid may "rebounce" , thus preventing black-hole formation 0. It has been argued, on the other 
hand, that a collapsing spherically symmetric charged dust will inevitably produce naked singularities due to shell- 
crossing The latter indicates the intersection of matter flow-lines along certain spacelike hypersurfaces. Although 
these singularities are considered weak [?|, since the curvature invariants and the tidal forces remain finite, their 
appearance could also signal that a nonzero Lorentz force and charged spherical collapse are physically incompatible. 
This tentative 'conjecture' is supported by the results of the present paper, which focuses on the role of the magnetic 
Lorentz force during the gravitational collapse of charged matter. 

We consider the gravitational contraction of an inhomogeneous spatially flat Tolman-Bondi model, filled with a 
pressureless highly conductive fluid, and allow for a perturbing weak magnetic field. This describes to a collapse 
of a weakly magnetized charged plasma. The weakness of the field is measured by its contribution to the total 
energy of the system. Our aim is to investigate the fate of the contracting fluid in the presence of a nonzero and 
without any symmetry constraints Lorentz force. To the best of our knowledge this question has not been addressed 
analytically. Nonzero Lorentz force means that, despite the absence of fluid pressure, the particle worldlines are 
no longer timelike geodesies. As a result, we expect the initial spherical symmetry of the collapse to break. We 
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quantify the consequences of the magnetic presence by looking at the convergence of neighboring particle worldlincs. 
This is done by analyzing the magnetic contribution to the Raychaudhuri equation. The latter describes the volume 
evolution of a given fluid element and has played a fundamental role in numerous studies of gravitational collapse and 
also in the formulation of the major singularity theorems (e.g. see ^3)- A relatively weak magnetic field contributes 
to the Raychaudhuri equation primarily via its Lorentz force. The benefit of our approach is that it identifies the 
effects of Lorentz force on the collapsing matter directly. In particular, the magnetic input splits up into a pair of 
stresses one of which always supports against the gravitational pull of the matter. Based on the weakness of the field, 
we ignore its backreaction on the Tolman-Bondi metric and use the latter to evaluate the aforementioned magnetic 
stresses. We find that the magnetic presence triggers as range of effects with an overall impact that depends on 
the specifics of the field in a rather involved way. These effects can severely distort the spherical symmetry and 
completely dominate the advanced stages of the Tolman-Bondi collapse despite the low levels of the magnetic energy 
input. Interestingly, of all the magnetic stresses, we found that those supporting against the collapse grow faster and 
we have also identified physically plausible magnetic configurations where this happens. Although one should be very 
cautious before extrapolating a linear result into the nonlinear regime, our analysis seems to agree with earlier work 
on magnetized collapse Put together, these studies suggest that a nonzero Lorentz force may be physically 

incompatible with spherically symmetric collapse and that there might exist situations where the support of the field 
can outbalance the total gravitational attraction, at least along certain directions. 



II. WORLDLINES OF MAGNETIZED MATTER 



Assume a general spacetime filled with a highly conductive perfect fluid and allow for a magnetic field. High 
conductivity means that there is no electric field and that the magnetic field is 'frozen in' with the matter. This is 
the well known MHD approximation (e.g. see 0). In what follows we will investigate the implications of the pure 
magnetic component of the Lorentz force on the collapse of such a model. We will do so by testing the convergence 
of the particle worldlines using the covariant approach to general relativity [Tl) . 

Covariantly, the dynamics of gravitational collapse is monitored through the Raychaudhuri equation, which describes 
the volume evolution of a self-gravitating fluid element. Consider a congruence of timelike worldlines tangent to the 
4-velocity field u a (with u a u a = —1) that follows the motion of the fluid. Raychaudhuri's formula determines the 
evolution of O = V ra it a , the scalar measuring the average contraction (or expansion) between a pair of neighboring 
particle worldlines [h| . In a magnetized environment we have |12| 

9 + ie 2 = -i (p + 3p + B 2 ) -2(<j 2 - oj 2 ) + D a u a + ii a u a , (1) 

where p and p are respectively the energy density and pressure of the fluid, B 2 = B a B a measures the energy density 
and the isotropic pressure of the magnetic field (B a ), a 2 and ui 2 are the respective magnitudes of the shear and the 
vorticity associated with u a and u a = u b V b u a is the 4-acceleration. When the right-hand side of the above is negative 
definite, an initially converging family of worldlines will focus (i.e. Q — > — oo) within a finite amount of time |10| . 
Thus, positive definite terms in the right-hand side of |QJ will resist against further gravitational contraction. 

The magnetic contribution to the Raychaudhuri equation comes form the energy density of the field which adds to 
the gravitational pull of the matter and also from the magnetic input to the 4-acceleration. The latter satisfies the 
momentum-density conservation law, which for a magnetized, highly conductive perfect fluid takes the form [T^ 

(p + p + §5 2 ) ii a = -D aP - e abc B b cm\B c - U ab ii b . (2) 

In the above D a = h a b V a is the covariant derivative operator orthogonal to u a , with h„h representing the projection 
tensor (i.e. h ab u b ), while Ii ab = —B( a B b } describes the anisotropic pressure of the field [l^- The second term in the 
right-hand side of Eq. is the Lorentz force, which is always normal to the field vector. Note that we consider 
non-geodesic worldlines, since the motion of the particles is dictated by the combined Einstein-Maxwell field and not 
by gravity alone. Also, the fluid flow is generally not hypersurface orthogonal, which explains the presence of the 
vorticity term in Q). 



III. TOLMAN-BONDI COLLAPSE 



The general inhomogeneous collapse of pressure-free spherically symmetric matter is monitored by means of the 
Tolman-Bondi metric, which in the synchronous gauge reads 

ds 2 = -dt 2 + R' 2 {1 + kr 2 )- 1 ^ 2 + R 2 dn 2 , (3) 
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where R — R(r, t), dfl 2 — dO 2 + sin 2 8d(j) 2 and a prime indicates differentiation with respect to r (e.g. see 01). The 
spatial curvature index k = 0, ±1 corresponds to flat, closed and open geometry respectively. Here we will only 
consider the k = case and throughout this paper we will use geometrized units with k — 1 — c. In the adopted 
coordinate system, the energy momentum tensor of a dust cloud is T a b — p(r, t)u a Ub, with u a = (1, 0, 0, 0) representing 
the 4-velocity of the fluid. Solved on the above metric Einstein's equations give 

R = R(r,t)=(~) M^ir-tf 3 and p = p(r, t) = —— , (4) 



AJ y ' 8ttR 2 R> 

with the functionals M — M(r) and r = r(r) describing the spatial distribution of the matter at fixed time. To 
proceed further we use the residual gauge freedom on the radial coordinate to define M — pr 3 /3, which makes p the 
mass per unit coordinate volume (i.e. 4irr 3 /3). Thus, the only physically free function left to describe the shape of the 
pressureless fluid at a fixed time is r(r). Note that the point t = r corresponds to a singularity which is reached at a 
different time for each shell of fixed radius r. Indeed, following and the Ricci scalar R a a = —8irT a a — 8np 
associated with any given shell of radius r = r* , will diverge when t = r(r*). Here, we assume that the arrow of time 
increases from t = 0, which marks the beginning of the collapse, to t = r(r). This in turn guarantees that r — t > 
for each r. Then, the density and the shear magnitude of the Tolman-Bondi solution are 

9 ~ 2tt(t - t) [3(r - t) + 2rr'] " ~ 3[3(t - t) + 2rr'] 2 (r - t) 2 ' U 

respectively. As we approach the singularity the behavior of these two variables changes and for t — > t their evolution 
is monitored by the following approximate expressions 
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- and a ~ , , (6) 

47rrr'(r-i) 9 (t - tf W 

indicating that the shear can dominate the final stages of the collapse. Note that, according to we guarantee a 
positive definite energy density for the matter by demanding that r' > at all times. The latter also offers a sufficient 

condition for avoiding cross-shell singularities (see 0] for further discussion). Finally, we stress out that by using the 
null coordinates, with 

dR = R'dr + Rdt and dv = dt - (R - l^dR , (7) 

the line element © reads 

ds 2 = -(1 - R 2 )dv 2 + 2dvdR + R 2 dn 2 . (8) 

Then, the surface where the radial velocities of null congruences vanish (i.e. for dR/dv = 0) is a spherically symmetric 
horizon. Therefore, by putting ds/dv = dR/dv — dfl/dv — 0, the horizon is the surface R 2 = 1. As we are looking for 
collapsing solutions only (i.e. with R < 0), the horizon will be on the restricted surface Rr = —1 or, more explicitly, 
on t(r H ) = T(r H ) - 2/^/9. 



IV. MAGNETIZED TOLMAN-BONDI COLLAPSE 



Observations have long established the widespread presence of astrophysical magnetic fields, while compact stellar 
objects are capable of supporting considerably strong fields. Neutron stars, for example, can carry magnetic fields 
that reach up to 10 G and 10 G. Despite their strength, however, the energy density of these B-fields is much 
smaller than that of the supporting matter (i.e. B 2 /p <C 1). One can therefore use its relative weakness to treat the 
magnetic field as a perturbation on a matter dominated background. Here, we will employ the Tolman-Bondi metric 
to study the collapse of an inhomogeneous, weakly magnetized dust cloud. We will do so by adopting the familiar 
MHD approximation, which is supported by the expected very high conductivity of stellar interiors. 

According to , the volume evolution of a weakly magnetized fluid element within the irrotational Tolman-Bondi 
spacetime is monitored by the following version of the Raychaudhuri equation 

G + ±<d 2 = -±p-2a 2 +D a u a + u a u a , (9) 

given that p = and B 2 <C p. The last pair of terms in the right-hand side of the above is entirely due to the 
magnetic presence, since ii a = e a b c B h cxa\B c / p for p = (see Eq. and also 0|). In what follows, we will employ 
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the Tolman-Bondi metric to evaluate these two terms, while ignoring the magnetic backreaction on the shear and the 
vorticity of the background model. 

At the MHD limit the electric field vanishes and Maxwell's equations reduce to a set of one propagation and one 
constraint equation. In covariant form, these are given by the respective expressions 

h a b £ u B b + OB a = and h a b V a B h = , (10) 

where £ u B a = u b \7bB a — B b VbU a is the Lie derivative of B a along the 4-velocity of the fluid. Adopting spherical 
polar coordinates and recalling that B a u a = 0, we set B a = (0, B r , B e , B^) with B a = B a (t, r, 6, (j)) and a = r, 9, <\>. 
Then, solving i|l(Jh) on our Tolman-Bondi background we obtain 

_ f_a (11) 

(r - t) [3(t - t) + 2rr'} ' V ; 

with F a representing time-independent functionals (i.e. &tF a = 0). Similarly, written relative to a spherical polar 
coordinate system, constraint IjlUb ) translates into 

sin 6(rF' r + 2F r ) + r(sm Qd B F B + cos 6F e ) = . (12) 

As with the matter density and the shear magnitude before, the behavior of the magnetic field changes as one gets 
closer to the singularity. In particular, for t — > r, we find that the magnitude of the above given i?-field evolves as 

1 AF 2 



where B 2 = B a B a and A — 3 -\/4^ 2 /81. This in turn combines with result (jHt.) to provide a measure of the energy- 
density ratio B 2 j p near the Tolman-Bondi singularity 

r2 AF 2 rr' 

T-V^W- (14) 

with A = it A. Since the magnetic density grows faster than that of the collapsing matter (compare l|13fl to Eq. ljU<)), 
the above also allows for a rough upper bound on the ratio B 2 j p. It is therefore clear that the weak- field approximation 
(i.e. B 2 j p <C 1) will hold as long as 

r-t » (AF?rT') 3/5 , (15) 

The r-dependence in the right-hand side of this condition implies that the weak-field approximation can be satisfied 
at any time during the contraction. For instance, for any arbitrary finite value of r — t, one can always set < r <C 
(t — t) 5 / 3 / AF 2 t' to ensure that (|15|) holds (recall that r' > 0). In other words, we are always able to identify a shell, 
labeled by its radius r, where condition l|15f) is satisfied. 

Remaining within the weak-field limit, we will now employ the Tolman-Bondi metric to evaluate the magnetic 
impact on the contracting spacetimc. Focusing on the later stages of the collapse (i.e. allowing t — ► r), the dominant 
components of last two terms in the right-hand side of the Raychaudhuri equation (see ©) read 

5 AF 2 t' . . a 4 AF 2 (dgF r ) 2 T 2 

D «■ - -3 and g 9 (r-*)"/3 ' (16) 

respectively. Thus, as long as F r ^ and deF r ^ the former of the above assists the contraction and the latter acts 
against it. Note that both terms grow faster as we approach the Tolman-Bondi singularity, compared to the matter 
density and the shear magnitude given by l|fik,) and (JHJd) respectively. In this case the global spherical symmetry 
of the collapse will be destroyed by the Lorentz force. Note that, as the supporting magnetic stress l|16b ) is the 
fastest growing near the singularity, the overall contribution of the field will tend to resist the collapse and this 
could cause the converging worldlines to bounce. In addition, when cross-shell singularities are not allowed, halting 
the gravitational contraction of a shell of proper radius r — ro means that all shells with r > ro will also cease 
collapsing. This is possible while still within the weak-field limit because condition l|15|l can be satisfied arbitrarily 
close to the singularity. In other words, although the magnetic energy density associated with a given collapsing 
shell has negligible contribution to the total energy-momentum tensor, the Lorentz force dictates the symmetries of 
the collapse. The latter results from the generic inhomogeneity of the Tolman-Bondi spacetime. It is still possible, 
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however, that nonlinear contributions to the metric, mainly caused by the non-geodesic motion of the matter, could 
overwhelm the magnetic resistance and push the system into further (anisotropic) contraction. 

In what follows we will demonstrate that fields with the aforementioned properties that support against gravitational 
collapse can be obtained as solutions of Maxwell's equations on the Tolman-Bondi background. Indeed, imposing the 
condition dgF r ^ on the S-field and then solving Eqs. l(TU|) we arrive at the functionals 

f f 
F r = --cos9 and F e = -sin0, (17) 

2 r 

where / is a constant. Therefore, both components of the associated magnetic field have a clear (^-dependence. Also, 
despite the fact that Bg diverges at r = 0, the energy density of the field is perfectly regular there because B 2 oc gggFg 
and ggg oc r 2 . Following (jT7Jl . however, the radial component of the field and its 0-derivative vanish at 9 — ir/2 and 
at 9 — 0, 7r respectively. Hence, along these directions the supporting magnetic effect disappears and the contraction 
will proceed uninhibited by the presence of the field. This behavior, which is probably typical (see |l6||l. could be 
seen as a direct consequence of the generically anisotropic nature of the field. The result is an extremely distorted 
collapse. In particular, while certain directions will continue collapsing, the gravitational contraction of most of the 
magnetized particles will face strong resistance by the Lorentz force. 

At this point we should also emphasize that the (9-dependence of the radial magnetic component is crucial for the 
future of the converging worldlines. To be precise, when dgF r = the last term in Eq. © approaches the expression 

1 AF 2 (F 2 + F 2 sin 2 9) 

u u a ~ ry 8 t L CIS) 

UaU ~ 4 ( T -i) 8 / 3 ' K ' 

instead of Qlfifa ). In this case the later stages of the magnetized collapse are dominated by Qlfifa ,) and therefore the 
contraction will proceed unimpeded. 

Finally, let us consider the homogeneous limit r' — > 0, which corresponds to FRW geometry. In this special case the 
shear contribution to the Raychaudhuri vanishes, while p oc (r — i)~ 2 (see Eqs. (0). Also, for r' = the right-hand 
side of both expressions in (|16[) vanish, which means that one needs to evaluate the magnetic stresses at higher order. 
Then, the total magnetic effect near the initial singularity is given by the sum 

u a u a + D a u a ^K{r-t)- % / 3 , (19) 

where K = K(r, 9, (f) is an involved function of the F a s and their derivatives and vanishes for an homogeneous 
magnetic field. The interesting point is that, as t — > r, the above stress grows faster than the matter density and 
therefore it is expected to dominate the final stages of the collapse. In addition the functional K does not have a 
definite sign and therefore its effect on the collapsing dust depends on the particular magnetic configuration. 



V. DISCUSSION 



Past studies of magnetized gravitational contraction have indicated that the presence of the field could affect 
the outcome of the collapse in nontrivial ways. In £|, for example, magnetism was found capable of halting the 
contraction of a finite cylinder, while in 5] the authors provided a static solution of the Einstein-Maxwell equations 
with a pressure- free matter component. It has also been shown that the spherically symmetric collapse of a charged 
star could bounce , although the bounce seems to cause naked shell-crossing singularities || . More recently, it was 
argued that magneto-curvature tension stresses could also affect the collapse of a conventional magnetized fluid 0] . 
All these claims suggest that the Maxwell field could play a key role during the gravitational collapse of a bounded 
system. Motivated by that we have considered the inhomogeneous contraction of a weakly magnetized Tolman-Bondi 
spacetime filled with a highly conductive pressureless fluid. The advantage of the Tolman-Bondi metric is that 
it offers the most general mathematical framework for studying dust collapse. Here, we did so by analyzing the 
magnetic contribution to the Raychaudhuri equation, which monitors the convergence of the particle worldlines in a 
covariant manner. We show that the magnetic input splits into two stresses, one of which always supports against 
contraction. Assuming that the energy density of the field is only a small fraction of the matter density, we have 
used the Tolman-Bondi metric to evaluate the input of the aforementioned stresses to Raychaudhuri's formula. Such 
an approximation was unavoidable given the fully analytic nature of our study. Within this limit, we found that 
the magnetic presence can severely distort the spherical symmetry of the collapse. Interestingly, we found that the 
magnetic stress which resists the collapse is the fastest growing one and we identified physically plausible magnetic 
configurations which allow this to happen. It should be noted, of course, that by ignoring the backreaction effects, 
mainly those due to the non-geodesic motion of the magnetized fluid, we have limited the range of our results. 
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Nevertheless, the tendency of the Lorentz force to dominate the collapse, even when the energy level of the field 
is relatively low, should not depend on the approximation level. This argument is supported by earlier studies, 
showing that spherically symmetric charged collapse produces naked singularities Thus, given that all known 
stars support magnetic fields of various strengths, we believe that these fields can play a protagonist's role in the 
evolution of such gravitationally bound systems. 
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